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We  describe  the  design  of  a  silicon-based  source  for  radiation  in  the  0.5-14  THz  regime.  This  new  class  of 
devices  will  permit  continuously  tunable,  milliwatt  scale,  cw,  room  temperature  operation,  a  substantial  ad¬ 
vance  over  currently  available  technologies.  Our  silicon  terahertz  generator  consists  of  a  silicon  waveguide  for 
near-infrared  radiation,  contained  within  a  metal  waveguide  for  terahertz  radiation.  A  nonlinear  polymer  clad¬ 
ding  permits  two  near-infrared  lasers  to  mix,  and  through  difference-frequency  generation  produces  terahertz 
output.  The  small  dimensions  of  the  design  greatly  increase  the  optical  fields,  enhancing  the  nonlinear  effect. 
The  design  can  also  be  used  to  detect  terahertz  radiation.  ©  2008  Optical  Society  of  America 
OCIS  codes:  040.0040,  040.6040,  060.4080,  130.0130,  130.2790,  130.3120. 


1.  INTRODUCTION 

The  applications  of  terahertz  (THz)  radiation,  roughly  de¬ 
fined  as  electromagnetic  radiation  with  a  frequency  rang¬ 
ing  0.5-10  THz,  are  a  topic  of  active  research  [1],  They  in¬ 
clude  medical  imaging  [2],  remote  sensing  [3],  chemical 
detection  [4] ,  and  tactical  imaging  [5] .  One  significant  im¬ 
pediment  in  this  field  is  the  expense  and  complexity  asso¬ 
ciated  with  THz  sources.  Current  methods  for  generating 
THz  frequencies  include  using  spectrometers  to  filter 
blackbody  radiation  [6] ,  pulsed  lasers  combined  with  non¬ 
linear  optical  materials  to  exploit  difference  frequency 
generation  [7,8],  quantum  cascade  lasers  [9],  and  high- 
frequency  microwave  mixers  [10]  among  others.  Another 
recent  approach  has  been  to  combine  quantum  cascade  la¬ 
sers  and  difference  frequency  generation  [11].  Some  re¬ 
cent  work  on  THz  devices  has  focused  on  metamaterials, 
which  may  enable  switching  and  control  of  THz  radiation 
[12],  but  this  has  not  yet  resulted  in  significant  progress 
on  new  THz  sources.  All  of  the  aforementioned  methods  of 
generating  THz  radiation  suffer  from  either  very  low  op¬ 
erating  temperatures,  lack  of  tunability  or  very  low  power 
output,  or  can  function  only  in  pulsed  mode.  The  majority 
of  current  approaches  to  THz  generation  in  fact  suffer 
from  several  of  these  problems  [13]. 

Current  approaches  to  using  y2  nonlinearities  suffer 
from  low  conversion  efficiencies  [7,8].  This  is  due  in  part 
to  the  fact  that  THz  modes  must  have  relatively  large  lat¬ 
eral  dimensions,  as  well  as  the  fact  that  nonlinear  fre¬ 
quency  conversion  is  inherently  less  efficient  when  the 
output  is  at  lower  frequencies.  That  is  why  current  ap¬ 


proaches  generally  require  pulsed  lasers.  Several  recent 
developments  have  made  possible  a  new  approach  to  dif¬ 
ference  frequency  generation,  however.  First,  it  has  be¬ 
come  possible  to  guide  optical  modes  with  relatively  low 
loss  in  small,  high-index-contrast,  ridge  waveguides.  Re¬ 
alizations  in  SOI,  for  example,  achieve  7  dB/cm  of  loss  in 
a  0.5 X 0.1  //m  waveguide  [14].  Second,  nonlinear  poly¬ 
mers  have  become  increasingly  powerful,  with  an  r33,  the 
standard  metric  of  nonlinearity,  of  170  pm/V  already 
demonstrated  [15].  It  is  expected  that  a  further  3X  im¬ 
provement  will  be  possible,  with  r33  values  that  will  ap¬ 
proach  500  pm/V  [16],  about  1  order  of  magnitude  stron¬ 
ger  than  typical  nonlinear  crystals  [17].  While  these 
nonlinear  polymers  are  typically  used  for  electrooptic 
modulation,  we  have  shown  unambiguously  that  the  non¬ 
linearity  is  ultrafast  [18].  When  combined  with  the  proper 
THz  waveguide  design,  and  with  the  relatively  large 
amounts  of  laser  power  that  are  now  available,  we  will 
show  that  practical  amounts  of  cw  THz  radiation  can  be 
generated. 

2.  THEORETICAL  BACKGROUND 

Maxwell’s  equations  can  be  written  for  a  nonlinear  me¬ 
dium  as  follows  [19]: 

/  s0dte  -  V  X  \  /  E  \  /-  E0dtxfjkEjEk 

\VX  p0dt  j\H/  \  0 
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Here  the  subscripts  i,  j,  and  k  represent  a  summation 
over  the  x,  y,  and  z  components  for  each  vector.  We  will 
consider  our  waveguide  to  be  described  by  a  dielectric 
function  that  is  continuous  in  the  propagation  direction, 
which  we  will  consider  to  be  the  2  direction.  The  inclusion 
of  metals  presents  no  difficulties  if  the  dielectric  constant 
is  permitted  to  approach  infinity  in  some  locations.  A 
guided  mode  will  always  have  the  form 


CiWif/iix^expiifaz-icoit).  (2) 


We  will  neglect  the  effects  of  loss  for  the  moment.  Physi¬ 
cally,  the  fields  are  always  real.  Therefore,  the  most  gen¬ 
eral  form  in  which  a  series  of  propagating  waves  in  a 
waveguide  can  be  written  is 


Hj  =  2  (Ci(z)ifii(x,y)exp(i/3iZ  -  ia-t) 

+  c*(2)i/f*(3c,y)exp(- 1/3,2  +  (3) 


Note  the  z  dependence  of  the  mode  constant.  For  a  linear 
material,  c  will  not  have  any  2  dependence.  But  for  a  non¬ 
linear  material,  this  will  not  necessarily  be  true.  Putting 
Eq.  (2)  in  Eq.  (1)  yields 


I  Hy(x,y ) 
-  Hx(x,y) 


2 


0 

-Ey(x,y) 

Ex(x,y) 

°  / 


exp(t/S/2  -  twit)  +  c.c. 


Xyjk 

xljk&j&k 

0 

0 

\  0 


(4) 


Here  the  mode  vector  has  been  expanded  to  show  the  in¬ 
dividual  components,  c.c.  denotes  the  complex  conjugate 
of  the  previous  expression.  Expanding  Eq.  (4)  yields  a  se¬ 
ries  of  different  equations  for  different  harmonic  compo¬ 
nents.  Let  us  make  the  discussion  more  specific.  Consider 
oq  and  olI‘2  to  be  two  optical  beams  near  1550  nm.  These 
beams  can  easily  be  selected  to  generate  a  difference  fre¬ 
quency  in  the  THz.  For  instance,  the  difference  frequency 
between  a  beam  at  1550  and  1558  nm  is  about  1  THz.  Let 
c d3  be  this  difference  frequency: 


4c3(z) 


/  \ 


3y 

-h3iX 

0 

-E, 


E 


3,  x 


exp(i/S32  -  ioi3t) 


\  °  / 

2e0(-  iio^expiifaz  -  i/32z)exp(-  ico3t)ci(z)c2*(z) 

ICex,kA 


XyjkEljE2,k 
X zjkElJE2k 

0 

0 

0 


(5) 


Here  f?1;  E2  and  so  forth  denote  the  relevant  electric  or 
magnetic  field  patterns  for  each  mode.  It  should  be  noted 
that  Eq.  (5)  is  not  quite  correct;  it  is  possible  that  addi¬ 
tional  modes  could  exist  at  co3,  which  in  principle  could  be 
excited  as  well  by  the  nonlinear  perturbation.  However, 
consider  that 


I  Hhy  \ 


I 


(E*x  E*  E*  H*  H*  H*  ) 


-Hi 


■/' 


Ej,X 

\  0  / 


dA 


(E  ■  x  H  +  E  x  H* )  ■  zdA . 


(6) 


It  can  be  shown  [20]  that  any  two  modes  at  a  given  fre¬ 
quency  with  differing  propagation  constants  are  orthogo¬ 
nal  under  the  metric  given  in  Eq.  (6).  Any  degenerate  sub¬ 
space  can  of  course  be  orthogonalized  with  the  usual 
methods.  As  a  result,  if  we  assume  that  there  is  no  degen¬ 
eracy  in  the  mode  set  involved  in  Eq.  (5),  we  can  integrate 
the  entire  expression  across  the  x,y  plane  to  solve  for  dzc3. 
It’s  worth  noting  that  if  the  modal  field  patterns  are  nor¬ 
malized  so  that  for  each  mode,  Eq.  (6)  evaluates  to  1,  then 
the  corresponding  physical  mode  as  found  in  Eq.  (4)  has  a 
time-averaged  power  flow  of  1  W.  We  adopt  this  conven¬ 
tion,  and  as  a  result  the  time-averaged  power  in  a  mode  in 
watts  is  equal  to  the  modulus  squared  of  the  correspond¬ 
ing  modal  amplitude  C;(z).  Integrating  Eq.  (5)  yields 

4c3(2)  =  -  2e0(-  iw^expd^z  -  i/32z  -  i/33z) 

Xci(2)c2(2)  J  E3  iy2ljkEl  iE2jdA.  (7) 


In  the  event  that  the  propagation  constant  difference  is 
zero  (/81-/S2-y33  =  0),  Eq.  (7)  can  be  written  in  a  simpler 
form.  Let  us  also  assume  that  the  only  important  nonlin¬ 
ear  contribution  in  Eq.  (7)  stems  from  the  nonlinear  mo¬ 
ment  along  a  single  lateral  direction,  taken  to  be  x  with- 
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out  loss  of  generality.  Finally,  the  possibility  of  loss  for 
mode  3  is  added.  In  this  case, 

dzc3{z)  =  -  2(-  iw3)ci(z)c2(z)x%xxy-  a?3c3(z),  (8) 


where  the  figure  of  merit  y  has  units  of  volt-seconds  per 
square  meter  and  is  defined  to  be 


y-£ 0  J  (9) 

Let  us  further  assume  that  the  two  input  beams  have 
negligible  loss  compared  to  the  generated  THz  beam,  and 
that  only  a  small  fraction  of  their  power  is  converted  to 
THz  radiation.  It  is  then  possible  to  write  down  the  final 
amount  of  power  that  is  generated  in  mode  3  due  to  non¬ 
linear  conversion  from  modes  1  and  2: 


P3(L)  : 


4co23PiP2(x2xxx?)2 
2 


-(1  -  exp(-  a3L))2.  (10) 


We  note  that  the  phase  matching  condition  in  Eq.  (7)  can 
be  rewritten  as 


n((Oi)Wi  —  7z(£02)Cl>2 

n3= - .  (11) 

£0 1  -  "2 

We  show  in  this  paper  that  perfect  phase  matching  is  at¬ 
tained  in  a  continuous,  nonsegmented  waveguide  by  engi¬ 
neering  the  three  effective  refractive  indices  of  the  com¬ 
posite  waveguide. 

Typically,  electro-optic  polymers  are  characterized  by 
their  r33  coefficient.  This  is  defined  as  follows  [17]: 

1  1 

“72 - 2=-®zr33-  (12) 

n7  n. 


That  is,  the  shift  in  refractive  index  along  the  poling  axis 
(taken  to  be  z  here)  is  defined  as  a  function  of  the  product 
of  r33  and  the  driving  electric  field.  This  can  be  repre¬ 
sented  as  a  x2  for  small  shifts  in  index  as 


v2 

Azzz 


n-ir33 

2 


(13) 


This  is  true  only  if  the  nonlinear  effect  is  not  attenuated 
at  higher  frequencies.  In  general,  this  is  often  untrue; 
however,  the  nonlinear  mechanism  involved  in  electro¬ 
optic  modulation  by  nonlinear  polymers  such  as  these  is 
based  on  a  deformation  of  an  electron  orbital  and  can  thus 
be  expected  to  be  ultrafast.  Additionally,  we  have  previ¬ 
ously  demonstrated  optical  rectification,  an  ultrafast  phe¬ 
nomenon,  based  on  these  polymers  [18],  and  our  results 
suggested  that  the  \2  will  be  relatively  unattenuated  at 
even  the  optical  frequency. 

The  refractive  index  of  these  nonlinear  polymers  is 
typically  around  1.7.  Therefore,  if  500  pm/V  of  tuning  is 
obtained,  we  can  expect  a  y2  of  2  X  10-9  m/V. 

Before  discussing  a  particular  waveguide  configuration, 
it  is  useful  to  calculate  the  order  of  magnitude  of  results 
that  might  be  expected.  Provided  that  a  given  interacting 
waveguide  geometry  achieves  the  phase-matching  condi¬ 
tion  [Eq.  (11)]  between  the  various  modes,  the  governing 
features  will  be  the  areas  that  the  modes  are  in,  their  ef¬ 


fective  indices,  and  the  extent  of  overlap  in  the  region  of 
the  nonlinear  material.  If  the  modes  are  normalized  so 
the  moment  in  Eq.  (6)  is  1  W,  then  the  average  electric 
field  will  be  approximately 


yp  /  Mo  \ 1/4 

\l2An  \  eo  / 


(14) 


Here  P  is  the  normalized  power,  which  we  set  to  1  W  by 
convention,  while  n  is  the  effective  index  of  the  mode. 
Given  this  expression,  one  can  then  calculate  Eq.  (9)  as 
follows: 


/  Mo 

I 

\  Eo 


3/4  1  Ac  1 

^ni?i2n3  yA1A2A32 


y p. 


(15) 


Here  A3,  A2,  and  A3  are  the  areas  of  the  relevant  modes; 
n i,  n2,  and  n3  are  the  effective  indices;  and  Ac  is  the  area 
where  all  of  the  modes  and  the  \2  material  overlap.  P  is 
the  normalized  time-averaged  power  of  all  three  modes, 
which  is  always  set  to  1  W,  included  here  for  clarity  of 
units.  With  this  expression,  it  is  easy  to  estimate  what 
conversion  efficiencies  might  be  approximately  obtained 
for  a  given  interacting  waveguide  geometry.  Let  us  imag¬ 
ine  that  we  have  engineered  a  situation  where  A3=A2 
=  10-13  m2,  A3  =  84  X  10-12  m2,  and  Ac=  1  X  10-13  m2,  and 
all  effective  indices  are  2.  In  this  case,  y  will  be  8.8 
X  10-4  V  s/m2. 

As  will  be  shown,  confinements  of  THz  modes  to  these 
kinds  of  area  are  possible,  though  the  waveguides  are 
lossy.  In  one  instance,  the  loss  near  1  THz  is  about 
25  dB/cm.  In  this  case,  ay  =  287  m_1.  Assuming  that  both 
optical  beams  are  at  a  power  level  of  0.5  W,  the  output 
power  would  be  1.5  mW.  This  level  of  output  power  would 
be  reached  to  within  1  dB  in  about  7.5  mm  of  device. 


3.  PROPOSED  DEVICE 

As  one  will  conclude  from  Section  2,  the  challenges  of  de¬ 
signing  a  silicon  terahertz  generator  break  down  into 
three  areas:  First,  designing  waveguides  that  will  guide 
both  optical  and  terahertz  radiation,  both  without  crip¬ 
pling  losses;  second,  optimizing  the  modal  overlap  with 
the  nonlinear  materials  to  provide  efficient  difference  fre¬ 
quency  generation;  and  third,  the  waveguides  designed  so 
that  phase-matching  conditions  between  the  THz  and 
near-IR  (NIR)  signals  are  met. 

It  is  well  known  that  silicon  can  be  used  to  guide  radia¬ 
tion  in  the  NIR,  and  that  silicon  is  relatively  transparent 
in  much  of  the  THz  [21].  In  addition,  because  of  silicon’s 
high  index,  silicon  waveguides  can  concentrate  optical 
modes  to  areas  of  1/1000  the  size  of  the  modes  of  ordinary 
optical  fibers  and  can  be  evanescently  coupled  to  nonlin¬ 
ear  optical  polymers  [22].  Silicon  microscale  waveguides 
can  be  efficiently  coupled  both  to  free  space  and  to  fibers 
[23],  and  large-scale  integrated  photonic  structures  can 
be  made  using  such  waveguides  [24,25] . 

Inexpensive,  high-power  sources  in  the  range  of 
1480- 1600  nm  have  recently  become  widely  available, 
owing  to  demand  from  telecommunications  applications 
[26].  Restricting  pump  sources  to  this  frequency  range, 
the  frequencies  available  as  a  difference  of  the  two 
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sources  range  from  dc  to  15  THz.  High-power  lasers  are 
also  available  in  specific  bands  in  the  1300-1480  nm 
range,  though  there  are  not  yet,  to  our  knowledge,  single 
sources  that  are  tunable  across  the  entire  range.  Based  on 
current  laser  and  erbium-doped  fiber  amplifier  technol¬ 
ogy,  one  can  expect  to  obtain  approximately  0.5  W  of  cw 
laser  power  [27] ,  though  this  power  can  be  expected  to  im¬ 
prove  in  coming  years. 

Our  proposed  STG  design  consists  of  a  silicon  core  for 
guiding  NIR  light,  within  a  metal  waveguide  for  confining 
THz  radiation;  optimal  performance  would  be  obtained 
from  copper,  but  other  metals  could  be  used.  This  basic 
geometry  is  shown  in  Fig.  1.  A  second-order  optical  non¬ 
linearity  enables  the  two  optical  signals  to  mix  and  gen¬ 
erate  a  difference  frequency,  and  various  dimensions  can 
be  optimized  in  order  to  achieve  phase  matching.  The  en¬ 
tire  structure  is  supported  by  a  micromachined  silicon  pil¬ 
lar. 

The  optical  waveguide  in  the  STG,  shown  in  Fig.  2,  con¬ 
sists  of  a  strip  of  silicon  within  a  dual-strip  metal  THz 
waveguide.  This  entire  structure  is  embedded  in  a  clad¬ 
ding  of  electrooptic  polymer  [28]  with  a  refractive  index  of 
1.7.  This  waveguide  geometry  creates  a  significant  eva¬ 
nescent  NIR  field,  which  overlaps  with  the  polymer  clad¬ 
ding.  Panel  A  in  Fig.  2  shows  the  waveguide  structure  and 
the  optical  mode  overlaid,  for  the  500-nm-wide  wave- 
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Fig.  1.  (Color  online)  Diagram  of  the  terahertz  generation  de¬ 
vice.  An  isometric  view  of  the  device  is  shown.  The  narrow  rect¬ 
angle  corresponds  to  the  silicon  waveguide  atop  an  oxide  pillar, 
while  the  base  corresponds  to  bulk  silicon  that  has  been  micro- 
machined.  The  thick  rectangles  indicate  the  metal  waveguide 
structure,  made  of  copper. 


A  0.5  um 
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Fig.  2.  (Color  online)  Diagram  of  the  optical  waveguide.  Dia¬ 
gram  is  shown  in  panel  A,  with  contours  of  the  mode  rendered  in 
increments  of  10%  of  \E\  assuming  1  W  of  input  power.  The  wave¬ 
guide  used  was  500  nm  wide  and  100  nm  tall;  the  width  of  this 
waveguide  is  varied  as  a  design  parameter.  The  Silicon  wire 
waveguide  is  supported  by  a  100  nm  wide  pillar  of  Si02.  In  panel 
B,  the  dispersion  diagrams  of  several  waveguide  widths  are 
plotted. 

guide.  The  effective  indices  of  refraction  from  three  differ¬ 
ent  Si  waveguide  widths  are  presented  in  panel  B.  The 
modes  were  solved  for  using  a  Hermetian  eigensolver 
[14].  With  such  guides,  phase  matching  for  THz  modes 
with  effective  indices  ranging  from  2.8  to  3.2  is  practical. 

Similar  optical  waveguides  have  been  shown  to  have 
losses  of  approximately  3.5  dB/cm  and  are  expected  to 
improve  significantly  in  the  future  [29].  From  perturba¬ 
tion  theory,  the  loss  for  the  optical  mode  from  interactions 
with  the  surrounding  metal  waveguide  can  be  bounded  to 
not  exceed  3  dB/cm,  and  it  is  generally  significantly  less. 
As  typical  device  lengths  for  the  nonlinear  interaction  re¬ 
gion  are  expected  to  be  below  2.5  mm,  optical  losses  of 
even  10  dB/cm  would  not  present  a  limitation. 

In  designing  the  THz  waveguide,  we  wished  to  achieve 
two  goals.  First,  we  wanted  to  make  the  THz  mode  more 
compact,  since  this  maximizes  the  figure  of  merit,  as  can 
easily  be  seen  from  Eq.  (15).  Second,  we  wanted  a  wave¬ 
guide  that  was  relatively  broadband.  Whereas  spanning 
as  much  as  10  THz  in  the  optical  regime  does  not  amount 
to  a  large  fractional  shift  in  frequency,  the  same  is  cer¬ 
tainly  not  true  in  the  THz  regime.  As  a  result,  we  sought 
a  design  inspired  by  a  transmission  line,  since  TEM 
modes  often  span  frequency  regimes  of  many  orders  of 
magnitude. 

To  span  the  maximum  frequency  range,  two  THz  wave¬ 
guide  designs  (labeled  I  and  II)  are  presented.  They  are 
both  similar  to  a  conventional  copper  dual-strip  transmis- 
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Table  1.  Description  of  Terahertz  Waveguides  and 
Device  Performances" 


Terahertz 

Waveguide 

B 

C 

D 

E 

F 

Type 

A  (/im) 

(um) 

(um) 

(um) 

(um) 

(um) 

I 

14 

3.5 

>20 

2 

2 

2 

II 

4.7 

3 

>20 

1 

1.3 

2 

^Columns  A  through  F  correspond  to  the  dimensions  in  Fig.  1 .  The  data  in  column 
C  refer  to  the  fact  that  this  dimension  is  to  exceed  the  specified  value.  Label  G  is  not 
present,  as  the  optical  waveguide  width  is  varied  independently. 


sion  line.  However,  the  portion  of  the  transmission  line 
where  the  mode  resides  is  partially  filled  with  a  plug  of 
silicon,  which  is  part  of  the  supporting  silicon  pillar.  The 
purpose  of  this  feature  is  to  raise  the  effective  refractive 
index  of  the  THz  mode  and  to  ensure  that  the  phase¬ 
matching  condition  is  met.  Precise  dimensions  of  the  two 
designs  are  given  in  Table  1. 

The  complex  geometry  of  the  THz  guide  makes  mode 
solving  challenging.  The  modes  were  solved  by  launching 
the  dual-strip  transmission  line  pattern  into  a  full  3D 
simulation,  and  observing  the  field  pattern  at  various 
points  down  the  waveguide.  A  steady-state  analysis  was 
used  on  the  Yee  grid  [30].  The  discretization  was  0.1  /xm, 
and  the  run  length  was  around  120  /mi.  From  the  relative 
uniformity  and  z  behavior  of  the  fields  so  obtained,  it  was 


clear  that  this  method  yielded  the  accurate  modes.  Once 
the  modal  pattern  is  known,  the  losses  and  effective  indi¬ 
ces  of  the  modes  were  easily  calculated  via  the  relevant 
closed-form  expressions.  Note  that  in  the  actual  simula¬ 
tion,  a  perfect  metal  was  assumed.  However,  using  per¬ 
turbation  theory,  the  waveguide  loss  due  to  imperfect 
metals  could  then  be  calculated. 

In  Fig.  3,  the  modal  patterns  for  THz  waveguides  I  and 
II  are  shown  at  various  frequency  extrema.  At  higher  fre¬ 
quencies,  the  THz  mode  tends  to  migrate  into  the  region 
with  more  silicon.  This  has  the  unfortunate  effect  of  pull¬ 
ing  the  field  away  from  the  optical  waveguide  and  de¬ 
creasing  the  generated  power.  The  effective  indices  of  the 
terahertz  waveguides  are  shown  in  panel  B  of  Fig.  4. 

The  waveguide  loss  of  the  THz  mode  is  highly  fre¬ 
quency  dependent.  This  is  due  to  the  effect  of  increased 
frequency  on  Ohmic  losses,  as  well  as  frequency  depen¬ 
dence  of  the  mode  pattern  The  oxide  pillar  supporting  the 
optical  waveguide  will  not  significantly  impact  the  THz 
mode  loss,  since  the  relevant  mode  overlap  is  very  small. 
Bulk  polymer  losses  as  high  as  30  dB/ cm  could  be  toler¬ 
ated  without  significantly  altering  device  performance.  It 
is  anticipated  that  future  nonlinear  polymers  will  easily 
meet  this  standard  [16].  Substrate  leakage  is  not  expected 
to  be  a  major  issue,  provided  that  the  support  pillars  are 
20  fin i  or  more  in  height.  Even  over  3-D  simulations  of 
100  fim  or  more  of  waveguide,  no  substrate  leakage  is  evi- 


Fig.  3.  (Color  online)  Terahertz  waveguide  modal  patterns.  The  modal  patterns  for  terahertz  waveguide  I  is  shown  in  panels  A  and  B, 
at  frequencies  1  and  6  THz,  respectively,  for  1  W  of  propagating  power.  The  modal  patterns  for  terahertz  waveguide  II  are  shown  in 
panels  C  and  D  at  6  and  14  THz,  respectively.  In  all  cases,  the  \E\  field  contours  are  plotted  in  increments  of  5%.  The  maximum  field  in 
Volts  per  meter  is  also  indicated  on  a  scale  bar.  This  is  for  a  mode  with  a  time-average  energy  of  1  W. 
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Fig.  4.  (Color  online)  Characteristics  of  terahertz  waveguide  modes  and  output  power.  Panel  A  shows  the  effective  index  of  the  terahertz 
modes  as  a  function  of  frequency.  Panel  B  shows  the  waveguide  loss  in  decibels  per  centimeter.  In  panel  C,  the  output  power  is  specified 
for  devices  specified  in  Table  1,  in  decibels  of  output  power  versus  the  output  frequency  in  terahertz. 


dent  in  such  a  configuration.  While  a  small  amount,  per¬ 
haps  less  than  0.1  dB,  might  remain  undetected  in  these 
models,  this  still  provides  an  upper  bound  of  10  dB/cm, 
which  is  far  less  than  the  losses  due  to  other  sources. 

Once  the  effective  indices  for  the  modes  are  estab¬ 
lished,  and  the  precise  modal  patterns  are  known,  it  is 
possible  to  calculate  the  conversion  efficiency  via  the  re¬ 
lations  presented  in  Section  2.  Table  2  and  Fig.  4  provide 
a  summary  of  the  cw  output  powers  achievable  as  a  func¬ 
tion  of  wavelength.  It  is  possible  to  change  the  waveguide 
width  lithographically,  suggesting  that  these  devices  may 
well  serve  in  an  array  to  maximize  the  frequency  coverage 


on  a  particular  chip.  Output  power  in  the  high  microwatt 
level  is  predicted  across  the  entire  frequency  range  from 
0.5  THz  up  to  15  THz. 

An  important  characteristic  of  these  devices  is  that  the 
generated  THz  beam  can  be  interchanged  with  one  of  the 
optical  beams.  This  can  be  seen  by  noting  that  Eqs.  (9) 
and  (10)  are  symmetric  in  the  electric  fields  of  all 
waveguides,  to  a  complex  conjugation  and  a  frequency 
factor  in  front.  As  a  result,  if  two  optical  beams  at 
0.5  watts  each  generate  a  1  THz  output  of  0.5  mW,  then 
an  optical  beam  at  0.5  W  and  a  THz  beam  at  0.5  mW  will 
generate  an  optical  beam  of  18  mW  in  the  exact  same  de- 


Table  2.  Full  Device  Designs  with  Their  Performance  Data0 


Device 

No. 

Optical 

Waveguide 

Width 

(rt 

Device 

Length 

(mm) 

Terahertz 

Waveguide 

Type 

Output 

Frequency 

Range 

(THz) 

Output 

Power 

Range 

(dBm) 

Input 

Laser 

Range 

Needed 

(nm) 

1 

0.375 

5 

I 

0.5-2 

-13.7,  -8.1 

1480- 

1600 

2 

0.4 

4 

I 

2-2.5 

-8.2,  -8.5 

1480- 

1600 

3 

0.425 

3 

I 

2.5-3 

-8.6,  -9.2 

1480- 

1600 

4 

0.45 

3 

I 

3-3.5 

-9.6,  -10.2 

1480- 

1600 

5 

0.5 

3 

I 

3.5-5 

-11,  -12.9 

1480- 

1600 

6 

0.6 

2.5 

I 

5-6 

-14.3,  -16 

1480- 

1600 

7 

0.4 

2 

II 

6-14 

-3.4-11.2 

1300- 

1600 

"The  two  input  lasers  are  assumed  to  be  500  mW  in  cw  power.  The  y2  value  is  taken  to  be  500  pm/ V,  or  2  X  10-9  m/V.  The  optical  waveguide  loss  is  taken  to  be  6  dB/cm, 
which  is  a  safe  upper  bound  on  expected  losses.  Note  that  the  optical  waveguide  width  corresponds  to  letter  G  in  Fig.  1 .  The  predicted  output  spectra  generated  by  these  devices 
are  plotted  in  Fig.  4. 
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vice.  That  is,  an  STG  can  be  used  to  shift  a  THz  signal 
into  the  optical  regime  and  also  add  about  15  dB  of  power. 
The  output  beam  will  be  located  at  the  pump  optical  fre¬ 
quency  plus  the  THz  frequency,  allowing  for  THz  detec¬ 
tion  schemes  that  will  use  conventional  diode-  or 
photomultiplier-tube-based  optical  detectors,  which  are 
extremely  sensitive. 

4.  HEATING  AND  FREE-CARRIER 
LIMITATIONS  ON  PERFORMANCE 

A  typical  problem  with  integrated  optical  CW  devices  is 
the  dissipation  of  waste  heat.  Clearly,  any  form  of  loss 
that  does  not  result  in  scattered  radiation  will  eventually 
convert  the  optical  input  energy  to  heat.  Under  continu¬ 
ous  operation,  this  heat  source  will  raise  the  temperature 
of  the  waveguide  and  polymer  until  it  is  balanced  by  heat 
dissipation  mechanisms.  If  the  steady-state  temperature 
exceeds  a  damage  threshold,  the  device  will  fail. 

A  thorough  discussion  of  the  extraction  of  waste  heat, 
and  the  limitations  this  will  impose  on  cw  input  laser 
power,  is  beyond  the  scope  of  this  paper.  Also,  such  con¬ 
siderations  will  be  driven  in  part  by  the  precise  character¬ 
istics  of  the  polymer,  parameters  that  can  be  engineered 
with  this  specific  function  in  mind.  However,  we  are  opti¬ 
mistic  that  the  power  levels  we  have  described  could  be 
supported.  We  have  already  demonstrated  that  cw  optical 
power  levels  of  100  mW  can  be  supported  in  such 
waveguides  clad  with  nonlinear  polymers  [22].  Typically, 
damage  thresholds  are  encountered  a  few  decibels  beyond 
this  level,  at  around  200  mW.  To  support  the  two  input  la¬ 
sers  of  0.5  W  that  we  have  proposed  above,  the  damage 
threshold  would  need  to  be  increased  by  less  than  1  order 
of  magnitude,  which  suggests  it  is  achievable  with  careful 
engineering. 

Another  consideration  is  the  fact  that  even  though  ra¬ 
diation  at  a  1300  nm  wavelength  and  longer  is  beyond  the 
silicon  band  edge,  various  mechanisms  can  produce  free 
carriers,  especially  at  high  power.  One  such  mechanism  is 
two-photon  absorption  (TPA)  [31].  In  very  small  silicon 
ridge  waveguides  such  as  the  ones  we  propose  to  use, 
surface-state  absorption  can  also  create  free  carriers  [32] . 
This  creates  nonlinear  loss  on  the  optical  signal,  which 
can  lead  to  rapid  signal  attenuation.  To  determine 
whether  this  would  be  a  limitation,  nonlinear  loss  mea¬ 
surements  were  performed  on  silicon-on-insulator 
waveguides  with  the  0.5  X  0.1  /xm  dimension  given  above. 
It  was  found  that  the  loss  within  a  waveguide  was  well 
modeled  by 

d I 

-  =  -a1I-  a2I2,  (16) 

d  z 

where  I  is  the  propagating  cw  intensity  of  a  beam,  and  at 
is  0.00014  yum-1,  while  a2  was  1.8  X  10-6  yum1  mw-1.  This 
corresponds  to  about  12  dB/ cm  of  additional  loss  at  a  cw 
power  of  150  mW.  While  this  is  certainly  problematic,  it 
should  be  remembered  that  the  THz  waveguide  losses  are 
far  higher,  and  so  the  interaction  length  in  which  a  THz 
beam  can  be  generated  is  fairly  short.  Therefore,  the  op¬ 
tical  losses  are  not  important  in  most  cases  until  they 
cause  significant  attenuation  in  only  a  few  millimeters. 
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All-optical  pump-probe  measurements  were  also  per¬ 
formed  in  order  to  determine  minority  carrier  lifetimes; 
this  work  will  also  be  detailed  in  a  forthcoming  publica¬ 
tion  [32].  The  minority  carrier  lifetime  was  found  to  be 
around  1  ns,  which  can  be  compared  to  the  typical  values 
of  1-10  ns  reported  by  Jalali  [33]. 

The  nonlinear  loss  can  be  reduced  by  removing  the  free 
carriers.  Several  methods  have  been  demonstrated  for 
electrically  sweeping  excess  carriers  from  a  silicon  wave¬ 
guide  to  reduce  effective  carrier  lifetime  [34,35].  Addition¬ 
ally,  it  has  been  demonstrated  that  ion  implantation  of 
Ar+  ions  can  reduce  the  tree-carrier  lifetime  to  100  ps 
[36],  which  would  amount  to  a  decrease  in  the  nonlinear 
loss  that  our  waveguides  experience  by  a  factor  of  10.  This 
approach  could  also  be  easily  integrated  into  our  process, 
since  it  does  not  require  any  change  in  the  waveguide  ge¬ 
ometry  or  electrical  contact.  Finally,  it  is  possible  to  treat 
silicon  waveguides  such  as  these  to  reduce  surface-state 
absorption  and  waveguide  loss  [37].  Lowering  the  wave¬ 
guide  loss  would  also  help  to  raise  the  thermal  damage 
threshold. 

Absorption  of  the  THz  mode  by  free  carrier  accumula¬ 
tion  is  expected  to  be  less  important.  This  is  due  to  the 
fact  that  though  the  absorption  coefficient  of  THz  radia¬ 
tion  due  to  free  carriers  is  approximately  100  times  larger 
than  the  absorption  coefficient  at  1550  nm  [38],  the  area 
of  the  silicon  waveguide  is  around  1%  of  the  area  of  the 
THz  mode,  suggesting  that  the  THz  mode  will  not  experi¬ 
ence  larger  losses  than  the  optical  mode  due  to  a  given 
free-carrier  concentration.  Of  course,  the  THz  mode  al¬ 
ready  experiences  losses  that  are  much  higher  than  the 
optical  signal  to  begin  with,  so  it  would  take  a  proportion¬ 
ally  larger  amount  of  loss  to  seriously  degrade  device  per¬ 
formance. 

While  free-carrier  losses  would  be  an  issue  in  using  the 
current  silicon  ridge  waveguides  in  a  device,  the  current 
nonlinear  loss  is  within  1  order  of  magnitude  of  what 
would  be  needed  for  the  proposed  designs.  Furthermore, 
there  are  numerous  examples  in  the  literature  of  solu¬ 
tions  to  this  problem.  It  is  therefore  likely  that  with  suf¬ 
ficient  engineering,  optically  induced  free-carrier  absorp¬ 
tion  does  not  pose  a  fundamental  limitation  to  the 
proposed  device. 


5.  CONCLUSIONS 

We  have  presented  a  design  for  a  device  that,  if  realized, 
would  allow  the  room  temperature  cw  generation  of  THz 
radiation.  Other  than  the  nonlinear  polymer,  our  system 
consists  of  materials  that  are  commonly  used  in  semicon¬ 
ductor  manufacturing  processes.  We  have  also  shown  that 
the  two  most  important  obstacles  to  high-power  perfor¬ 
mance,  thermal  damage  and  optically  free-carrier  excita¬ 
tion,  are  likely  to  be  surmountable  with  a  reasonable  level 
of  engineering.  Our  STG  can  function  at  room  tempera¬ 
ture  and  can  provide  a  level  of  output  power  that  is  far 
larger  than  what  comparable  processes  can  currently 
demonstrate.  We  believe  that  our  device,  once  realized, 
has  the  potential  to  make  the  THz  regime  vastly  more  ac¬ 
cessible  to  science  and  industry. 
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